function skyplot2(az,el,ns)
%SKYPLOT        Generate satellite azimuth/elevation angle plot


satn=1;pos=1;% pos->文字注记位置，pos=1->代表注记内容在卫星第一次出现的位置
% 遍历各卫星(按列存储各卫星高度角、方位角信息) GPS->[1-32],BDS->[96-158]

    ind=find(abs(az(:))>0);%非空信息列ind
    if ~isempty(ind)
        % [row,col]=size(ind);
        rho = rad2deg(el(ind));% 极径=90度-高度角(deg)
        theta =az(ind); % 极角=方位角(rad)
        polarplot(deg2rad(theta), deg2rad(rho), '*','MarkerSize',10,'Color',[0.6,0.5,0.2]);
        % pos=fix(median([1:row]));% 文字标注在轨迹中间
    % for i = 1:ns
	% 	leg{satn}=num2str(i,'C%02d');satn=satn+1;
    %     text(deg2rad(theta(i))+0.1, deg2rad(rho(i))+0.3, num2str(i,'%02d'), 'HorizontalAlignment',...
	% 	'center', 'VerticalAlignment', 'middle','FontName','Times New Roman');
    %     hold on
    % end
    end

rlim([0,90]); % 设置极径范围（0到90度）
thetaticks(0:30:360)
ax=gca;
ax.ThetaZeroLocation = 'top';
ax.RTick=[0,15,30,45,60,75,90];
ax.RTickLabel={'','15^{o}','30^{o}','45^{o}','60^{o}','75^{o}',''};
ax.RDir='reverse';
ax.ThetaTickLabel={'North','30^{o}','60^{o}','East','120^{o}','150^{o}','South','210^{o}','240^{o}','West','300^{o}','330^{o}'};
% legend(leg);
set(gca,'FontName','Times New Roman','ThetaDir','clockwise','GridLineStyle','--','Fontsize',12);
